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Recent results presented in the structure functions working group are briefly sum- 
marized for the following topics: The theoretical treatment of heavy quarks in 
structure functions, higher-order corrections for the leading-twist evolution (inclu- 
ding small-rc resummations), the present status of the proton's parton densities, 
and the impact of higher twists raj determinations of the strong coupling constant. 
The reader is referred to Part I til for accounts of the transition to the photopro- 
duction region, pion and photon structure results, and high-Q 2 phenomena. 

1 Heavy quarks in structure functions 

The charm structure functions, especially Ffi, have attracted considerable in- 
terest over the past years. Unlike in the fixed-target regime, F% makes up a 
sizeable fraction, up to about a quarter, of the total F^ in the HERA small- a; 
region. Despite being suppressed, it contributes significantly as well to the 
scaling violations in the kinematic range covered by the NMC data. 

At low scales, Q 2 « m 2 , F% is uniquely calculated from the light parton 
densities via the j*g — > cc Bethe-Heitler process and its 0(c? s ) corrections! 2 ), 
without invoking the concept of a charm parton distribution (we neglect here a 
possible intrinsic charm component, which seems to be relevant only at high a;). 
For Q 2 3> m 2 large logarithms appear in the coefficient functions C2.L, which 
may require a resummation. At x < 10 -2 these logarithms dominate C2 
already for Q 2 > 20 GeV 2 , but C'l only above 10 3 GeV 2 u. Previous leading- 
order results □ for this resummation have been extended to higher orders in 
ref. This leads to a high-Q 2 description in terms of four massless flavours, 
with the charm distribution uniquely specified by the light parton densities. 

The problem of the transition between both approaches, i.e., the construc- 
tion of a variable flavour-number scheme (VFNS) has also been addressed trQ. 
There seems to be agreement .that a unique construction does not exist, and 
indeed the prescription of refs.Lm differ at non-asymptotic values of Q 2 . Note, 
however, that this ambiguity concerns only the coefficient functions, as the 
parton evolution can be kept strictly massless without any loss of generalityQ. 
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Results for the VFNS prescription of ref.B are compared to the non-resummed 
calculation (usually called fixed flavour- number scheme, FFNS) in Fig. 1. The 
differences are typically 10% at small x, i.e., they are of the same size as the 
factorization-scale dependence of the FFNS calculationsB, shown for one value 
of Q 2 in Fig. 1 as well. Hence also the latter approach seems to be applicable 
at the present level of accuracy in the HERA small- a; region. On the other 
hand, somewhat larger effects are possible for other VFNS prescriptions U. 
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Figure 1: Left part: comparison of F£ as obtained in NLO from FFNS and VFNS 
calculationsB Right part: the factorization scale dependence of the FFNS predictions at 
one representative Q 2 . The light parton densities in both parts are not the same. 



Until now the extraction of from D* production at HERA has been per- 
formed only in the non-resummed approach (FFNS), as the required exclusive 
cross sections have been calculated only in this framework so farEl. However, 
first results have been presented from a VFNS event-genemtor Monte-Carlo 
program for semi-inclusive heavy-quark production in DIS lj. 

As is obvious from Fig. 1, the stability of the calculations deteriorates 
towards larger x, i.e., towards the threshold region. This effect is even more 
pronounced at lower Q 2 , and renders present calculations in the range of fixed- 
target experiments rather unreliable. A Sudakov resummation of leading and 
next-to-leading threshold logarithms, as discussed in ref. O, may lead to a 
stable framework also in this regime. 
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2 Higher-order corrections for structure function evolution 



We confine ourselves to the unpolarized singlet case, where new results have 
been obtained during the last year. The corresponding 2-loop (NLO) splitting 
functions are fully (at all x) known for about two decades. On the other hand, 
at O(a^) only the four lowest even- integer moments were determined so farl!3. 
A first step towards the full 3-loop anomalous dimensions, which are required to 
match the accuracy of present and forthcoming high-precision data, has been 
taken recently by calculating the finite terms of the 2-loop operator matrix 
elementsliia. Moreover a first partial result has been derived by means of the 
large- Nf expansion, namely the (a s /47r) 3 NjCG contribution to -P 3g E-3: 



P g 3 g loop (x,N]C A ) = - ^|87<5(l-a;) + (304 + 172a; + 208a: 2 ) In x 

32 

- 48(1 + x) In 2 x + 32 - — + 192(1 + x)[ip'(l) - Li 2 (x)] 

[l-o:] + 

+ -^——^(52 + 19a; + 52a: 2 ) ln(l - x) + ^ - ^ (236 + 47a; + 236a: 2 ; 
x 3x 

The first moments of this expression agree with the results of ref. E3. 

An alternative approach for the small- a; region has been to resum the most 
singular small- a; terms of to all orders in a s . For the gluonic splitting func- 
tions P gq and P gg these terms read c^*-(l/x)cts ln(l/x) k ~ 1 , and the coefficients 
c^ x were determined long ago as well. More recently also the leading contribu- 
tions to P qq and P qg were derived 113, which contain one power of ln(l/a;) less 
than their gluonic counterparts. These terms dominate the respective split- 
ting functions at some very low values of x, depending on the size of the less 
singular contributions. Until recently estimates of the impact of such terms, 
which tends to be enhanced substantially by the ubiquitous Mellin convolu- 
tion, were only possible by educated guesses based on momentum conservation 
constraints and the structure of the LO and NLO splitting functions t3. 

An inclusion of the leading small- a; logarithms into the analysis can lead 
to very good fits of all small- a; Fi and F% HERA data, as demonstrated in 
a scheme-independent evolution.-equation approach £3 all well as within the 
framework of k± factorization 113. However, both approaches seem to require 
values for observables and parameters - a rather small Fl in the first case and, 
more notably, a very small a s in the second one - which may be interpreted 
as phenomenological indications of large subleading corrections. 

During the past year the-calculation of 0(a s ) corrections to the BFKL 
kernel has been completed li*3'E2l. This result fixes the next-to- leading small- a; 
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Figure 2: The cumulative effect of the 
splitting function P gg (x,a s ) in the DIS s 
the first fifteen Lx and NLx iV-space exp 
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(NLx) piece of P gg . In the DIS scheme the presently known terms readEH 

OO , v I 

P g f S (N, a s ) = a s P gg , (N) + a 2 s P™?(N) + £ ( £ J + * , 

1=3 ^ / 

with a s — 3a s /7r and N the usual Mellin variable shifted by one unit. The 
Lx and NLx coefficients b ggi i are compared in Fig. 2, where also the resulting 
splitting function is shown for a s = 0.2. The NLcc corrections turn out to be 
exceedingly large in the HERA x-region, leading to a grossly negative splitting 
function already above x ~ 10~ 3 . Thus the ln(l/x) expansion is inapplicable 
to P gg at any x- values of practical interest, a situation which is by no means 
a special feature of the DIS scheme as demonstrated in ref . E3. 

Likewise, taking the NLx corrections at face value, the hard pomeron in- 
tercept lop would read wp(a s ) = 2.65 a s (l — 6.36 a s ) for Nf = 4E3, leading 
to negative values for Q 2 as high as about 300 GeV 2 . Hence a reliable exten- 
sion of structure-function evolution calculations beyond the usual fixed-order 
perturbation theory, if possible at all, will require new theoretical concepts. 
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3 Status of parton density parametrizations 

Major updates have been presented by the MRsli and GRV0 groups, super- 
seding their respective '96 and '94 parton setsEi The CTEQ collaboration 
has released a dedicated study of the uncertainty of the gluon density c3. The 
present (central) sets refer to values of a s (M|) = 0.116, 0.1175 and 0.114 for 
the CTEQ 4 Ell, MRST1 and GRV (98) distributions, respectively. Note also 
that heavy quarks are treated differently in these parametrizations. 

The present status of the gluon density is shown in Fig. 3. In the smalLa; 
part the preliminary HI and ZEUS error bands from F2 staling violations c3, 
as well as the recent HI results from DIS charm productionE3, are compared to 
these parametrizations. The difference between CTEQ4M ('massless charm') 
and CTEQ4F3 El (' massive charm') indicates the impact of the heavy quark 
treatment. The small- a; gluon density seems rather well constrained down to 
x ~ 10~ 4 . Note, however, that relevant theoretical uncertainties (estimated, 
e.g., by factorization-scale variations) are not taken into account in Fig. 3. 
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Figure 3: Present small-a: and large- x constraints on the proton's gluon distribution. 



Large uncertainties on xg still persist in the large- a; region, x > 0.2. Heie 
the classic constraint has been prompt-photon production in pp collisions l3. 
However, extractions of the gluon density from these data suffer from sizeable 
scale uncertainties as shown by the gray error bandE3. In addition, there is the 
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possibility of a sizeable gluon kx , as recently indicated by r esults from E706E1 
In fact, the MRST large- a; error bandE3 on xg stems from varying kx between 
(upper curve, set 2) and 640 MeV (lowest curve, set 3) in the fit to the WA70 
data. In view of these problems the CTEQ gluon-uncertainty analysis derives 
its error band from DIS and Drell-Yan data alone. It is interesting to note that 
all three bands are similarly wide for x > 0.3. By propagation to high scales, 
benchmark uncertainties have been derived for gluon-gluon and gluon-quark 
luminosities at the Tevatron and the LHC E3. 

As the total quark density is quite well constrained by F2 data for 10 -4 < x 
< 0.7, the other critical issues are the flavour decomposition and the x — > 1 
behaviour. In both areas new results have been reported. E&£6 has pub- 
lished their high-mass data on the pp/pd Drell-Yan asymmetry u. As shown 
in Fig. 4, these results strongly constrain the ratio d/u, especially in the range 
0.1 < x < 0.3. Preliminary results on this ratio, inferred from semi-irpjusive 
DIS, as well as data on F$/F% have also been presented by HERMES B. 
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Figure 4: The E866 Drell-Yan asymmetry data and their impact on the light-quark sea. 

F2/F2 is a dominant source of information on the d/u ratio. As FJ 1 is 
inferred from deuteron measurements, possible nuclear binding effects enter 
here Lj. Such effects have been neglected in recent parton parametrization. 
They can, however, lead to drastic modifications of the valence-quark ratio 
d v /u v at large x, as quantified in ref.E3. The region of very large x, x > 0.8, is 
usually not taken into account either. The DIS and resonance-region data of 
SLAC have been employed to study the absolute x — > 1 behaviour of the valence 
quarks. A significant flat contribution, as previously discussed in connection 
with the HERA high Q 2 events Ea, is found to be strongly disfavouredl23. 
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4 Higher-twist effects and the strong coupling constant 



Present-day determinations of a s from DIS structure functions involve rela- 
tively low values of Q 2 . Hence higher- twist corrections can have appreciable 
effects. A new result on a s (M 2 ) fimn the Gross-LLewellyn-Smith (GLS) sum 
rule has been reported by CCFREj. Their iron data, together with results 
from other neutrino experiments, in the region 1.3 < Q 2 /GeV 2 < 5.0 lead to 



NNLO 



(M Z 2 ) GLS - 0.114 



-0.009 



(exp.) ±0.05 (th.) 



The theoretical error is dominated by the uncertainty of the higher-twist c<j 
tribution. Nuclear 1 /Q 2 corrections have been studied for this case in ref. I 
Very small corrections to the GLS sum rule are found, but sizeable effects for 
the incomplete GLS integral, Sgls^t^O; Q 2 )i an d for the iron-to-nucleon ratio 
Rz{x,Q 2 ) = F 3 Fe /F 3 N , see Fig. 5. The latter results may be relevant for a s 
determinations from scaling violations in neutrino- nucleus DIScl. 
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Figure 5: Predictions for the nuclear effects in F£ I— 1 Left part: the incomplete GLS 
integral at Q 2 = 5 GeV 2 for A = d (dashed) and A = Fe (solid). Right part: the ratio 
R 3 for iron at Q 2 = 3 (solid), 5 (dashed), and 15 GeV 2 (dashed-dotted). 

New fits to the electromagnetic DIS data of BCDMS, SLAC and NMC at 
x > 0.3 have been presented as wellc3, including target-mass corrections and 
simple approximations for dynamical higher- twist terms. If the shape of the 
latter is allowed to vary freely, the fits result in 

a™ h0 (M 2 ) dQF *.«, = 0.114 ±0.002 (exp.) . 

Not surprisingly, this result agrees with previous analyses El. On the other 
hand, a very low value of a s {M 2 ) = 0.103 ± 0.002 (exp.) is obtained H3 if 
the sharjes of the twist-4 and twist-6 terms are adopted from the renormalon 
model Cj. Note that this finding seems to be at variance with the analysis of 
ref. E3, where good agreement between the large- a; data and the renormalon 
approach has been found for a fixed value of a s (M 2 ) = 0.120. 
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